function l = labornew(s)
% specify labor supply as a function of worker index s
% in the initial equilibrium calibrated to the US in 1980, labor supply is
% normalized to 1 for all s
% use this function for comparative statics with respect to labor supply

l = ones(size(s));

for i = 1:numel(s)
    if s(i) < 0.2
        l(i) = 0.9;
    elseif s(i) >= 0.2 && s(i) <= 0.4
        l(i) = 0.9 + (s(i)-0.2)*0.1/0.2;
    else
        l(i) = 1;
    end
end

%l = 1 - 0.2*(1-s).^4;
    
end